A low temperature derivation of spin-spin exchange in Kondo lattice model 
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Using Hubbard-Stratonovich transformation and drone-fermion representations for spin-| and for 
spin- 1 , which is presented for the first time , we make a path-integral formulation of the Kondo 
lattice model. In the case of weak coupling and low temperature, the functional integral over 
conduction fermions can be approximated to the quadratic order and this gives the well-known 
RKKY interaction. In the case of strong coupling, the same quadratic approximation leads to an 
effective local spin-spin interaction linear in hopping energy t. 
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I. INTRODUCTION 

In systems in which there is little direct overlapping 
of wave functions of localized magnetic atoms, the 
spin-spin interaction mediated by itinerant electrons 
is the dominant magnetic interaction, and has been 
an interesting topic for several decades. In particular, 
the study of spin-spin interaction in diluted magnetic 
semiconductors(DMS) has drawn much attention in the 
past several years. One of the key issues in explaining 
the magnetic properties of DMS is how the very diluted 
holes mediate the interactions between the localized, 
randomly distributed I -spins When the coupling 

between itinerant fermions and spins is weak compared 
to the Fermi energy, perturbation theory applies. Ruder- 
man and Kittcl initially considered the indirect exchange 
coupling of nuclear magnetic moments by conduction 
electrons[6|. And later, this interaction, conventionally 
called RKKY exchange, was extended by Bloembergen 
and Rowland^ , KasuyaJ^ and YosidaQ- Usually, 
RKKY interaction can be obtained using second order 
perturbation theory |Toj . But in the strong coupling 
regime, the conventional perturbation theory doesn't 
work and we need to devise a new way to derive the 
effective spin-spin exchange. In section II, we present 
a path-integral formulation of Kondo lattice model, in 
which local spins are represented by drone-fermions. 
Here the drone-fermion representation of spin-| is 
presented for the first time. At low temperatures, the 
spin-spin exchange is obtained. In section III we discuss 
a two site system and the origin of ferromagnetic inter- 
action linear in t. The last section is a brief summary. 



II. THE PATH-INTEGRAL AND THE 
SPIN-SPIN EXCHANGE 

We consider Kondo lattice model with Hamiltonian 



where cl a (ci a ) denotes the creation(annihilation) opera- 
tor of the itinerant electrons (or holes) and n\ a = c^cv 
is the number operator. Si denotes the localized spins, 
either 5 = | or 5 = §. s = he* tc is the electron spin 
where r = (ti , T2, T3) is the vector of the three usual Pauli 
matrices. Hopping energy ty = t > if i,j are nearest 
neighbors and zero otherwise, Jk > is the Kondo cou- 
pling and (i denotes chemical potential. The symbol (i, j) 
implies that the summation in the first term is taken over 
nearest neighbors. The other two summations are taken 
throughout the lattice. Note that there exists a local spin 
on every site and these local spins do not interact with 
each other directly. To derive the mediated interaction, 
we need to separate the conduction fermions and the lo- 
cal spins first. One way to this is to introduce an auxil- 
iary field and use Hubbard-Stratonovich transformation. 
For this, one needs to express the coupling S • s in terms 
of squares of Hcrmitian, bounded operators. There are 
different ways to do this. For example, one can write 

S-s = i[(S + s) 2 -S 2 -s 2 ] = -i[(S-s) 2 -S 2 -s 2 ] (2) 

Using S 2 = 5(5+1) and s 2 = -f (n T +n{) 2 + §(ra T +n;) 

S-s = I[(S+s) 2 -5(5+l) + J(n T +n i ) 2 -^(n T +n i )] (3) 
or 

S.s = -i[(S-s) 2 -5(5+l) + J(n T +n i ) 2 -^(n T +n i )] 

To formulate the path-integral, we also need a repre- 
sentation for the local spins. For spini,we can use the 
drone-fermi representation^: S z = p f — \, S + = 
fi(d + eft), 5~ = (el + eft)f where /, el are fermionic 
operators. One of the most important conveniences of 
this representation is that unlike the Schwinger fermion 
representation, this representation does not require con- 
straints. For spin- 1 , we have a generalization (to our 
best knowledge, this generalization is new, ref.|l2( only 
represents the spin operators in terms of fermi opera- 
tors and the representation is obviously not a drone- 
fermion representation since one term is bosonic and the 
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other is fermionic) : S+ = V3f±(d + S) + 2f\f 1: S~ = 
VS{d+S)h+2flf 2l S z = /J/i +2/J/2-I . Note that we 
use 3 fermi fields for spin-|. In general, this type of fermi 
representation is viable for S— spin when 2"" 1 < 2S+ 1 < 
2™, as counted in^l- Thus it is seen that for both spin- 
i and spin— | , we can write the partition function as 



purely fermionic path-integral. Using expression (3) and 
the relation (the so-called Hubbard-Stratonovich trans- 
formation) e~ A2 = -j= f^^dxe"^ 2 ^ 12 ^, which holds 
for Hermitian and bounded operator A , we can write 
the partition function as 



Z = j MiM^iWe-Z^S,!^^) J Vciai r)Vcl(r) J Vd^Vdl^Vf^Vf^r) 
where T> denotes functional measure and Ai(r), ifi(r) are auxiliary random fields and the extended Hamiltonian 

^ g /g 

H = t c i<r c J CT ~ l J K)^2(nii +n u ) +^[i v / 2JK(S i + s ; ) ■ A; + iJ - Jk^Kt + n n)} (6) 

(i,j>,<7 i i 

I 

Using the Fourier transformation A m k = bosonic fields, the electron part of the functional-integral 
N~ 1 Y,iP~ 1 ' 2 Jo ' dre^e^Aiir) where uj® = is written as 
2^7r, (m = 0,1,..) is the Matsubara frequency for 



no- 
where u n = 2 "^ l " 1 7r, (n — 0, 1, ..) is the Matsubara frequency for fermionic fields. Using the path-integral technique for 
Grassmann variable, we have 



lnZ c [A,(^] = In — — - — 

Dek[uj n d mn ) 

3 

= lnDct[(<9 r +tk — [i— -Ji<)<5p,k] — lnDet(9 T ) 



3/3 / t/x 

-lnDet[<5 mn (5 ki p + (iu m + t k - \i - -Jr) -1 ^ ^y^KV'm-n.k-p + iJ ^^ cr ' An-n,k- P )] (8) 



At low temperatures such that tJJk/P <C the largest of t, Jk and /i, the third term can be calculated perturbatively. 
Using InDetA = TrlnA, we have 



3/3 / </k 

Tr[(iw m +tk - M - 4 J K)~ 1 (jy ^ J KiPm- n ,k-p + iW ^Fg " ' An-n.k-p)] 



i Tr[(iw m + < k - A* - |^k) 



3 / </k 

;JKfa-n,k-p + ^'W dato ' 7 ' -^rrt-w.k-p)] +••■ (9) 



2Nf3 r ' v V 2N/3 
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To the quadratic order, i.e., neglecting the part • • • , we have 



lnDct[c5„ m <5 k ,p + {iu m + tu - \i - \ J K ) 1 (h/ 2 tv/3 2 yTO ~" ,k ~ p + y 27V/? 2 °" ' ^-"^-p) 



2W„ + £ k - ,U - |J K 
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1 ■ ^ . ^ *y' 2JV0 JviVm— ri,k— p + *y 2Nj3 <J ' ^-m— n.k— p iy 2N fi ^K'/'ri— m,p— k + *y 2Nfj (T ' ^n— m,p— k 
" P^k ^ ^ + *k - M - f^K «Wn + tp - M - | Jk ^ 

The first term vanishes in the thermodynamic limit N — > oo. Therefore the leading term is the second term 

x -» x 2N/3 ^K^m— n,k— p^Pn— m,p— k ~l~ 2N f3 ' 7i,k— p ' -^-n— m,p— k 

(^m + ik-M-f^KXiWn+ip-M-f^K) 

= — 5Z(3v? m ,ky-m,-k + ^m,k ' ^-m-k)Qmk (11) 
m,k 



where 



Qmk (iw m+ „ + t k+p -/i- |j K )(iw„ + ip-/i- |J K ) ^ ^ 



so the partition function 
2i 



x cxp{-^ ^(3^ ro , k ^, k + A m;k • < ik )g mk }e-^S mt s rak .A^, k (13) 



m.k 



where GL4 mk <i(/?„ k are just ordinary integration measures as A mk , <^„ k are now Fourier components. We have used 
replacement A mk — ► i^4 m k,^* nk — * ->^mk) i-c, redefine the Fourier component such that 

^k=^E^y /3 ^ T ^w 

Hence the induced spin-spin interaction comes from the integral over A, which gives 

J PA mk e- E ^ Amk - A ->-^ E ™^ A ™- k - A ™.^™ k -^S m Js mk .A^, k +A m , k .s^ k ) = e ^ Y: ™- kS ™ k i + 4Q mk Smk 

(14) 

Now since 

Q ii ^ [ ^ ^ ] (i^) 

„ * P - ik+p - i^m „ «w m+Il + t k+p - /i - |J K iu) n + t p -iJ,- |J K 



I 

, using the Matsubara frequency sum y^*L_„ ■ 1 , = where f (x) = d J z i is the Fermi function. In the low 
e _j^ +1 , we have temperature limit, we can take lu^ ~ so Q mk ~ Q k , 

n _£v- /(tk+p_ ^ _ t J K)-/(t p -M-!^K) 



^ „ tp - t k+p - iw^ 



p 



(16) 



4 



where 



Qv 



N ^ 



- V - p K ) - fjtp - y - %Jk) 



(In the case ik = k"/(2m), this is nothing but the func- 
tion -^3^)0 and we reproduce the RKKY interaction). 
If t S> Jk, the induced interaction is 



drVF(i-j)S r Sj 



J 2 

K \ "* c * 



mk 



Qk'^mk (18) 



m,k 



which is RKKY-like since it is proportional to J K and 
inversely to t which corresponds to Fermi energy. 

In the strong coupling case as in diluted magnetic semi- 
conductors, t < Jr. Keeping only the quadratic term of 
the functional determinant in (10), the resulting spin- 
spin effective action at low temperatures is of the form 



^E^k^ + Qk)- 1 ^^- J (M-t^JySi-Sj 
m,k s ij 



if Qk 3> j^- for all k , where Jij is the Fourier trans- 
form of the (J£ + Qk) _1 - This interaction is linear in t. 
The linearity in t of the induced spin-spin interaction is 
intuitive since when Jk is very large, then Jk is irrele- 
vant and the only energy scale left is t. A more formal 
discussion of induced spin-spin interactions for a two site 
sysytem is given in the next section. 



III. PERTURB ATIVE DERIVATION OF 
MOMENT-MOMENT EXCHANGE 

(1) Exchange between classical moments at half-filling 



We now consider a half-filled lattice with a clas- 
sical moment Sn on each site, where n is a unit 
vector. The relation of this system to the Kondo lattice 
model described by Hamiltonian (1) is as follows. The 
quantum-mechanical states of the local spin S can be ex- 
pressed in terms of the eigenstate |n) of S, S|n) = Snln). 
We can then obtain the energy levels by considering 
different orientations of n. In the strong-coupling limit, 
the hopping term is taken as a perturbation while the 
on-site interaction is the unperturbed Hamiltonian. 
Using the canonical transformation 



f 



s (m-ai + (n 1 «-l)A) (20) 



V2 - 2n iz 

*l = /0 ((l-nujod+m+A) (21) 
V2 - 2n iz 



where n\± = n- lx ± in\ y , and ai,/3i are new fermi fields, 
the on-site term becomes 



(19) JKS^ Si (r). ni (T) = lj K sY,[aHT)a i (T)-f3Ur)P i (r) 



(22) 

Hence at half-filling (the number of particles equals the 
number of lattice sites) , the unperturbed ground-state 
is (for Jk > 0) 



|Gnd) =I]>i|0} 



(23) 



with energy -Ecnd = —^JkSN. N is the number of lat- 
tice sites. Since the hopping part reads H' = 
where 



iJij = frj- — {[m-n i+ + (1 - n iz ){l - n Sz )]a\a s + [n ; _(n jz - 1) + (1 - nu)nj_]a[/3j 

2 y/(l-n iz ){l-n iz ) 

+ [{n iz - l)n,+ + n i+ (l - n iz )](3\a- } + [{n iz - \){n- JZ - I) + m+n^^} (24) 
the only non- vanishing term in Jf'|Gnd) is, 

Y^tiA 7= \. A n-fa* - 1) + (1 - n^nj-JaJ/JjlGnd) (25) 

Defining states 

|ij) = |/3l,.../3 i t ^,...,( j ),...,/3j v ) (26) 

in which site j is empty, site i is doubly occupied and all other sites are singly occupied by (3 quasi-particles. Therefore 
the second-order perturbation is 

(gndj giT^W E(ij) gjjlGnd) 

kl 
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Note that 



(ij|ff'|Gnd) = (-lf+H ii: 



- n iz )(l - n- iz ) 



[«i-(nj 2 - 1) + (1 - ni z )n- } 



(28) 



we get the second-order correction to the ground-state energy 



[ni+{n- )z - 1) + (1 - n iz )n 



j+J 



(29) 



AE, 



(2) 
Gnd 



tt, 



4J K S 3 3 Z^9Jt^.V v 



'■^2J K S y S 2 



(30) 



We see that this is an antiferromagnetic interaction, which agrees with references |14) and For 
half-filling, the ground-state is highly degenerate and we must use degenerate perturbation theory which 
not practical since the degeneracy is often too high to be handled. Eq.(3) indicates that the effective 
exchange which is linear in t obtained in section (II) applies only to systems away from half-filling. 

(2) Two sites energy-levels 

In the case of two sites, the Hamiltonian (1) is 



less than 
is in fact 
spin-spin 



/ Jk qz Jk q- 

/ _2 1 2_ 1 



H 



with eigenvalues 



For J K S > 2t, they are 



,t J J J 



Jk c 
t 



V o 




E 2 
E 3 
Ei 



±-\Jit 2 + J^S 2 ± 2J K SW2 + 2ni • n 2 



*f - 1^/2 + 2m- n 2 - - 1(2 + 2m ■ n 2 )] 



JkS , |i| 



^V2 + 2m n 2 - — [1 - -(2 + 2m • n 2 )] 



M_j£l V 2 + 2n 1 .n 2 + -|^[l-i(2 + 2n 1 .n 2 ) 



JkS , |t| 



2 V2 + 2ni ' n2 ' J K 5 



[l-i(2 + 2n 1 -n 2 )] 



(31) 



(32) 

(33) 
(34) 
(35) 
(36) 



If there is only one electron, only E\ is filled in the 
ground-state with energy 

Eo = -M_j£l V 2 + 2n 1 .n 2 --^[l-i(2 + 2n 1 .n 2 )] 

(37) 

and the correction due to hopping is 



1*1 



AE ~ -nV2 + 2m n 2 



(38) 



which is always ferromagnetic. When there are two elec- 
trons, i.e., the system is half-filled, Ei,E% are filled in 



the ground-state and the energy is 



t 1 



E = E 1 +E 2 =-J K S-—(l-n 1 -n 2 ) (39) 



So the correction due to hopping is 

2J K S" 



AEn 



: i-„ 1 .„ 2 )=^^L(„ i .„ j -l) (40) 



This agrees with the conclusion (30) for a general half- 
filled lattice. It is seen from (38) and (40) that whether 
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the induced moment-moment interaction is linear or 
quadratic in t depends on the number of particles. This 
justifies our previous argument. 



IV. CONCLUSIONS 

In this paper, we present a path-integral formulation 
for Kondo lattice model with localized spins [S = | or 
S = |). In both cases, the system can be described 
in terms of fermions conveniently since spin operators 
are difficult to deal with in either path-integral or Feyn- 
man diagram calculations. In the low-temperature limit, 
the conduction fermions can be integrated perturbatively. 
and the resulting effective spin-spin exchange is RKKY- 
like interaction in the weak-coupling regime. In the 
strong-coupling regime, the effective spin-spin exchange 
is linear in t if Qk ^ 7^ for all k where Qk is the function 



-p3(k)0] in the lattice case. Since Q(k) depends on fx 
which is the chemical potential, we see qualitatively that 
whether or not the leading order of the induced spin-spin 
interaction is linear in t depends on the particle density 
in the system. At half-filling, the leading term of inter- 
action is antifcrromagnetic and quadratic in t according 
to quantum-mechanical perturbation theory. It should 
be noted that though the path-integral formalism agrees 
with quantum-mechanical perturbative analysis qualita- 
tively, the former does not show that half-filling is of 
criticality. Therefore further study on the path-integral 
formalism is still called for. 
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